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Abstract 

In this paper we continue the study of the truncated conformal space approach 
to perturbed boundary conformal field theories. This approach to perturbation 
theory suffers from a renormalisation of the coupling constant and a multiplica- 
tive renormalisation of the Hamiltonian. We show how these two effects can be 
predicted by both physical and mathematical arguments and prove that they are 
correct to leading order for all states in the TCSA system. We check these results 
using the TCSA applied to the tri-critical Ising model and the Yang-Lee model. 
We also study the TCSA of an irrelevant (non-renormalisable) perturbation and 
find that, while the convergence of the coupling constant and energy scales are 
problematic, the renormalised and rescaled spectrum remain a very good fit to the 
exact result, and we find a numerical relationship between the IR and UV cou- 
plings describing a particular flow. Finally we study the large coupling behaviour 
of TCSA and show that it accurately encompasses several different fixed points. 
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1 Introduction 



The Truncated Conformal Space approach (TCSA) of Yurov and Zamolodchikov ||] is a 
widely-used method to study the finite-size dependence of perturbed two-dimensional con- 
formal field theories. It is based on truncating the infinite dimensional Hilbert space to 
a finite-dimensional system on which the Hamiltonian is studied numerically.[]lt has been 
known for a long time that the method has various convergence problems which can reduce 
its effectiveness j3j. The principal problems are a renormalisation of the coupling constant, 
a renormalisation of the energy scale and differences between ground state contributions in 
different sectors, all of which depend on the size of the truncated system (we shall always cut 
off the size of system by taking all states whose unperturbed energy above the ground state is 
less than or equal to a given number, which we shall call the truncation level). In a previous 
paper we showed how the renormalisation of the coupling constant in perturbed boundary 
conformal field theories could be studied using a variant of standard perturbed-conformal field 
theory methods [^]. In this paper we extend this study to the second effect, and show how 
the leading energy scale renormalisation is an overall multiplicative renormalisation which 
can be considered equivalent to a renormalisation of the size of the system. We find 'physical' 
arguments based on the operator product expansions and also more rigorous arguments based 
on an analysis of the eigenvalues of the perturbed Hamiltonian and show that these give iden- 
tical results. We test our results using two integrable conformal field theories, the tri-critical 
Ising model and the Yang-Lee model, as in both cases the finite-size spectrum has been found 
using TBA methods and this provides an accurate quantitative check of the proposed results. 

In the next section we briefly review the case of the tri-critical model to illustrate the 
renormalisation issues to be solved and to find numerical estimates for the coupling renormal- 
isation and energy rescaling. We then show that these have scaling forms and find numerical 
estimates for the associated exponents. 

In the third section we show how these can be derived, to one loop order, from considera- 
tions of the renormalisation of the perturbed action and in the fourth section show how these 
can be proven, to leading order, for all energy levels, from analysis of the eigenvalues of the 
perturbed Hamiltonian. In the subsequent two sections we check these predictions against 
the numerical data in the tri-critical Ising model and the Yang-Lee model 

In section || we consider the case of an irrelevant (non-renormalisable) perturbation and 
in section ^ we consider the flows beyond the fixed points and speculate on the exponents 
that have been found numerically. Finally in section [H] we present our conclusions. 



2 The TCSA approach and its errors 
2.1 The TCSA approach 

We start with a CFT defined on a strip < y < L of width L in the upper half plane with 
coordinate z = x + iy. We take the strip to have conformally invariant boundary conditions 
so that the system is conformally invariant, the Hilbert space carries an action of the Virasoro 
algebra and decomposes into a direct sum of representations of the Virasoro algebra. The 
representations occurring depend on the boundary conditions on the strip. 

1 In this paper we study the original form due to Yurov and Zamolodchikov, not the revised version of 0, H|. 
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The unperturbed CFT Hamiltonian generating translations along the strip is 

H= [ L T XX ^. (2.1) 
Jo 27r 

We will map the strip to the upper half plane with coordinate w = exp^nz/L) in terms of 
which the CFT Hamiltonian is 

H = i( L °-^)< < 2 - 2 > 

where Lq is the zero mode of the Virasoro algebra. 

We are interested in perturbations by one or more boundary fields <fii(x) living on the 
bottom edge of the strip, y = 0. We take these to be quasi-primary fields of conformal 
dimension hi. If the coupling to these fields are \ii then the perturbation is given by an 
addition to the action 

SS = I ^2ia<f>i(x)6x . (2.3) 

•* i 

When mapped to the upper half plane this gives the perturbed Hamiltonian as 

(2.4) 



H= ~L 



where yi = 1 — hi. Note that (|2.4j ) is only correct if the fields <j>i are primary; there are correc- 
tions if they are quasi-primary but not primary. We will normally consider the dimensionless 
operator 

« = (;)* = fa -a) m 

i 

where A, = H%L y are dimensionless coupling constants. 

The TCSA approach is to restrict this Hamiltonian to a finite dimensional space of ex- 
citation level n or lower. If we take the projector onto this space to be P n then the TCSA 
Hamiltonians are 

Hn = ( L ° " 2l) + £ ^ P nUm ■ (2.6) 

i 

These Hamiltonians can be diagonalised numerically and their eigenvalues and eigenstates 
form the TCSA approximations to the perturbed system. Any quantity that can be considered 
in the perturbed system can also be considered in the TCSA system; the only question is how 
good the approximation is, and whether the dependence on the truncation level is either small 
or can be estimated efficiently. In the next section we present the two leading truncation effects 
- the renormalisation of the coupling constant and the rescaling of the energy levels — in 
the case of the tri-critical Ising model. 



2.2 The tri-critical Ising model numerical results 

We take as our example model example the boundary tri-critical Ising model on a strip. 
We shall review the boundary conditions and their flows in section [j| for the moment it is 
sufficient to know that, amongst others, there are conformal boundary conditions labelled 
(11), (21) and (12), and the (12) boundary condition can be perturbed by a field of weight 
3/5 with the following flows: 

(11) ^_ (12) _ ► (12) . (2.7) 
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The spectrum of the Hamiltonian on a strip with an unperturbed (11) and a perturbed 
(12) boundary condition has also been analysed from the TBA approach by Feverati and 
collaborators in J?j. Using the results in plj] and [22, p3| , we can relate the TBA and conformal 
perturbation theory parameters and so we can compare the TCSA spectrum to the exact 
perturbed conformal field theory spectrum for this system. 

The most immediate difference is that there is an overall shift in the energy levels. For 
finite perturbations for which there are no divergences in the perturbation expansion this 
corresponds to a measurable free energy per unit length. In this case the perturbation theory 
is divergent and the overall shift in the energies is not meaningful: in TCSA it depends on 
the truncation level and does not converge to a finite value. As a consequence, in figure [j] we 
just show the difference between the ground state and the excited states - the energy gaps: 
in figure 1(a) we give the gaps as calculated using the TBA method, and in figure 1(b) the 
TCSA gaps. 




Figure 1: The energy gaps for the strip with boundary conditions (11) and (12) + A0( 13 ) as 
given by the TBA and TCSA methods plotted against A. 

On inspection of figure [l], one can also discern the two other effects. In the TCSA plot, 
there is an apparent convergence of many levels at or near A = 2. This corresponds to a fixed 
point of the flow where the (12) boundary condition has reached its (21) fixed point. This 
is only reached at A = oo in the TBA flow. The two couplings are related by a non-trivial 
renormalisation. Less obvious but still measurable is the overall rescaling of the energy levels. 
At A = —2 the TCSA spectrum has almost regrouped into equally-spaced levels but the 
spacing is approximately 15% larger than the (correct) spacing in the TBA spectrum. This 
effect grows markedly with increasing |A|. 

We can find numerical estimates for the renormalisation (of the coupling constant) and 
rescaling (of the energy scale) by comparing the TCSA and TBA data. Firstly we consider 
only the energy gaps to remove the unphysical ground state energy contributions in the TCSA 
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scheme. Next, since both the overall scale and the value of the coupling are different, we 
consider the ratio of two gaps to find a quantity which is independent of the energy rescaling 
and use this to determine the effective coupling constant. We can then use this in turn to 
find the energy rescaling. 

If we denote the coupling constant in TCSA by X n and that in the TBA by X^, and the 
n th energy gaps by A™ and AJ^ respectively, then we can calculate 



Ale A^ 
/oo(Aoo) — .1 ) /n(A n ) = -ry 



(2- 



We can then find the effective (TBA) coupling corresponding to a given TCSA coupling from 
the function g n defined by 



-^oo — foo{fn(X n )) — X n g n {X n ) 

and the energy rescaling r n (A n ) can then be found from 



AJo(Aoo) 



(2.9) 



(2.10) 



In figure ^| we plot these functions for the flows starting from the (12) boundary condition for 
three different truncation levels. It is the aim of this paper to show how these two effects can 
predicted from physical and mathematical arguments. 



2.5 - 





(a) The renormalised coupling constant X n g„(X n ) 
plotted against X n 



(b) The energy rescaling function r n (A n ) plotted 
against \ n 



Figure 2: The numerical coupling constant renormalisation and energy rescaling found for 
the tri-critical Ising model at truncation levels 6 (red, solid), 14 (green, dotted) and 22 (blue, 
dashed). 
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3 Physical Arguments 



3.1 The TCSA renormalisation group 



In |12j] , we presented an argument deriving the leading term in the coupling constant renormal- 
isation based on standard considerations of the partition function and the operator product 
expansion of the perturbing field. We repeat this argument here in a little more generality, 
both as a derivation of the coupling constant renormalisation which we shall check in the 
next section, but also to show that this does not in fact give the leading term to the energy 
rescaling, as might have been thought was the case. 

We consider a strip geometry of width L with complex coordinate z = x + iy with < 
y < L. We take a general boundary perturbation of the form 

^2^i4>i(x) , (3.1) 

i 

where /U, are the couplings to the boundary fields (pi which have conformal weights hi- We 
define y% = 1 — hi. 

We consider a truncation to level n which we can treat using a projection operator P n 
onto states of excitation level n or below. We can use the projection operator to restrict 
the perturbation to the truncated Hilbert space so that the TCSA approximation to the 
perturbation |3.1| is then given by 

Y,Vi(Pn<t>iPn) ■ (3.2) 

i 

The effect of the perturbation is to introduce the following operator 

Pexpl-/ ^2fn(P n (f>i(x)P n )dx) , (3.3) 

into all expectation values, where V is path ordering. Mapping the strip to the upper half 
plane by w = exp(7rz/L), this operator becomes 

Pexp ( - [°° J2HPnU^Pn)j^- ) , (3.4) 
y J x =o ^ (arc-)* J 

where Aj = ii%L Vi is the dimensionless coupling constant. 
Expanding this exponential out to second order we get 

x roc ^ x 

^ J x =o {xK)y* 

+ ^XjXk / P n ^{x)P n ^)P n \y . . , , yh + 0(A 3 ) ■ (3.5) 
^ J x=0 Jx'=o (xtt)* (x'vr)^ 

We define the TCSA renormalisation group by the requirement that the partition function 
is unchanged when the truncation level is increased. This introduces a level dependence into 
the coupling constants, which we now denote by Aj(n). 
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Due to translation invariance, we do not need to compare the whole expression in equation 
|3.5| , instead we can just compare the integrand with respect to x, at x = 1, which is 

£ *£U*(Dft "Eff £„ ^(DflAMft^J + O(A') . (3.6) 

We can further simplify matters by assuming that the Hilbert space contains a vacuum state 
|0) and a state \(f>i) corresponding to the field which satisfy 

(0i|^(l)|O) = Sij , <<^(l)0 fc (x)|O> = (1 - x) x ^- x *C ijk . (3.7) 

Sandwiching equation |3.6| between the states (4n\ and |0), we get 



^i( n ) - Aj(n)A fc (n) / 

jk J *=° 



;i)P n ^(x)|0)(7rx)^-^-^dx + O(A d ) . (3.8) 



This must be invariant under changes in n so that two second order in A. Using the fact that 
(0fc|0j(l)(P n +i — P n )(f)j(x)\0} is the coefficient of x n+1 in the expansion of (1 — x ) hl ~ h i~ hk and 
is equal to 

T{h 3+ h k -h l+ n+l) „»■ •> - 

Y(h j + h k -h i )T{n + 2) r(hj + h k -hi) { +U{ /n}) ' {6 - y> 
we find the TCSA renormalisation group equations 

~ n(Al(n + 1} " Ai(n)) = ? ^(^"^ r(fe .+f _^. } + °( A3 ) • ( 3 - 10 ) 

jk 3 l> 

3.2 The coupling constant renormalisation 

In the case of the perturbation by a single field with self-coupling C, the TCSA renormalisation 
group equation fl3.10| ) becomes 

HA C 

n-r- = -tt^ . A 2 + (9(A 3 ) , (3.11) 

with solution 

A(oo) = ^ . (3.12) 

1 " yr(h)(nn)y A(n) 

There are two predictions from this calculation - firstly the general prediction that the coupling 

renormalisation has a particular scaling form, namely that the function g n in equation 2Si 
has the form 

9n{x) = g{xn~ y ) , (3.13) 
and secondly a particular prediction for the 1-loop behaviour, 

9w =( i -sf^r + ° (i2 '""° ) <3 - i4) 



We test these predictions for the tri-critical Ising model in section |5.1| . We will find that 
the functions log(g n (— xn y )) are indeed almost identical for various values of n and in good 
agreement with the prediction, confirming both the prediction of the scaling form and the 
approximate numerical expression for this scaling function. 
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3.3 The energy rescaling 

We can try to apply the same arguments to derive the energy rescaling, assuming that this 
too comes from the operator product of the perturbing fields giving a contribution to the 
bare action (corresponding to the kinetic term in a standard Lagrangian theory). The Bare 
Hamiltonian is (tt/L)(Lq — c/24). The important part is the term Lq which is a mode of the 
energy-momentum tensor T(z). This field appears in the operator product of the boundary 
perturbing field as 

T(z)<P(w) ~ 2k <j>( w ) + -J— </>'(„) + ... , (3.16) 

[z — w) z z — w 

where we have also included the OPE of T with the perturbing field. 

If we denote the coupling to the perturbing field by A and to the energy momentum tensor 
on the boundary by A^ we then have = 2, y? = — 1 and we find the TCSA renormalisation 
group equations become 

dA T 2h/c 2 . / q1 m 

n— — = - — -z , —7 -X + . . . . (3.18) 

dn (n7r) 1 + 2 yr(-2y) v ; 

To find the leading dependence on the induced coupling At we can take A to be constant and 
using Ay(oo) = 0, find 

\ T (n) = - / — — dn 



dn 

2h/c A 2 



+ 0(A 3 ). (3.19) 



r(-2y)(l + 2y) (nvr) 2 ^ +1 

The addition of the term 

fj, T T(x)dx , (3.20) 



to the action will give the following term to the Hamiltonian (after mapping from the strip 
to the upper half plane) 

5H = u T T(l) = % (-—. \„ -(l - c/24) + other modes^i , (3.21) 

P KJ L\ T(-2y)(l + 2y) (rm) 2 y n\ 1 ) J 

so that the Hamiltonian becomes 

H + SH = 1 ( 1 - - 2h l° (l -c/24)+... , (3.22) 

L\ r(-2y)(l + 2y) {mrfv n J V 1 J ' v ; 

This gives the rescaling function as 

2h/c A 2 1 

rn(A) " 1 " T(-2y)(l+2y) {n^n + ' ' ' ' (3 ' 23) 



8 



Unfortunately, this has the wrong ra-dependence - it is one order of n too small. If we plot 
log(r n (A) — 1) against log(n) for fixed A for the functions shown in figure ||(b), in the tri-critical 
Ising model, we find that the prediction for the leading exponent 

log(r n (A)-l) = alog(n) + ... , (3.24) 

is approximately a = —0.85, which is close to — 2y = —0.8 and far from —1 — 2y = —1.8. As 
a consequence, we deduce that the coupling of the perturbing fields to the energy-momentum 
tensor is not the leading source of the correction to the energy rescaling function. 

Instead we look to the constant ground-state energy contribution coming from the coupling 
to the identity operator. As we show in the next section, this is actually not a constant because 
of the presence of the projector P n and this can also give a correction proportional to the 
bare Hamiltonian. 



3.4 Corrections from the identity operator 

In perturbation theory where the perturbing field has a weight greater than 1 /2 one normally 
ignores the ground state energy as it is a divergent unphysical quantity. In this case the 
TCSA cutoff makes the contribution finite, and further more, the presence of the projector 
P n means it is not constant. 

We are interested in the correction to the coupling to the identity operator which arises 
when the expression [O] acts on a state of excitation level E. If we denote this state by \E) 
and the coupling to the identity by \x(E, n), then we are interested in the expression 

(E\<p{l)P n (f){x)\E)—— (3.25) 



J x , =0 ' ' " '(x')y 
Requiring that this be invariant as we change the truncation level n, we get 
dAi(J5,n 



dn 



~ (Ai(£,n + 1)-Ai(£,n)) 
A(n) 



^t?) 2 f' 1 At 

/ W(l)(Pn+l - P n )m\E)-- , (3.26) 



where we only take the contribution to the OPE on the right hand side that comes from the 
identity channel. The projector (P n +i — Pn) picks out the state at level n + 1 in the action 
of the field <j)(x) on the state \E). The action of the field 4>{x) on the state \E) has a mode 
expansion 

^(x)\E)=Y,4> m x m - h \E) , (3.27) 

so that the term at excitation level n is just the coefficient of x n+1 ~ E . From the operator 
product expansion ( 0.15 ), the coefficient of x n+l ~ E in (E\<f)(l)(j)(x)\E} is 

T(n + l-E + 2h) 

r(n + 2 - E)T(2h) ' 1 ' 

This means that the leading contribution to the RG equation (|3.2(: for the coupling Ai is 
dAi(£,n) _ A(n) 2 f 1 T(n + 1 - E + 2h) x n+i-E dx 



dn TT 2 y J x=0 r(n + 2 - E)T(2h) (x)v 
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A 2 (n-E) 
T(2h) 
A 2 1 ' 



2/1-2 



+ ... 



E 



. 1 + 2y— + ...) (3.29) 
(nir) 2 y T(2h) \ n J y ' 

The first term is the constant contribution which we shall now ignore; the second is the 
excitation- level dependant term which we are interested in. We can integrate equation ( 3.29 ) 
using the boundary condition Xi(E, 00) = 0, to find 

»oo \ 2 



X 1 (E,n)\ E 



1 



A 



ir 2 v T(2h) 
1 



2yEn~ 2y - 1 dn 



(mr) 2 v T(2h) 



E 



(3.30) 



We can now replace the excitation level E by the operator (Lq — c/24), as E is approximately 
the eigenvalue of (Lo — c /24), and we finally find we have generated a term proportional to 
the bare Hamiltonian. We have two equivalent ways of understanding this term which lead 
to the same expression for the energy rescaling r n (A) 

Firstly, we can see that this coupling will lead directly to a change in the Hamiltonian; 
the new Hamiltonian is 

(1 + /.^LtvL XA) - c/24) + ±-6{l) 



7T 
I 



(mr) 2 f T(2h) 
giving the energy rescaling function as 



r n (X) = 1 + 



A 2 



1 



r(x) = 1 + 



+ 



(3.31) 



(3.32) 



(nvrfy F(2h) v ' TT 2 yT{2h) 

This has the correct n dependence (in agreement with the numerical data) and is also a very 
good fit to the actual rescaling function as we see in section 5.2 . 

An alternative viewpoint of the energy rescaling function is that it represent an effective 
change in the geometry of the system. Coordinate transformations are implemented in CFT 
by changes to the action (see e.g. [13 j). The coordinate change — > + corresponds to 
the change in the action 



5S 



1 

2^ 



T^a v d 2 x 



(3.33) 



The energy-dependant correction to the identity operator we have just calculated can also 
be put in this form. On the upper-half-plane, with complex coordinate w = rexp(iO), the 
correction is 



SS = -- [ Ai(£,n)— : 
vr J r 

Firstly, we note that on the upper half plane 



A- 



(nvr) 2 2T(2/i) 



dr 
Lc\ — . 



Ln 



(w 2 T(w) + w 2 T(w)) — 
=0 27r 



(3.34) 



(3.35) 



Combining equations (3.34) and ( 3.35| ) and transforming to the strip with coordinate z 
x + iy, the correction to the action becomes 

X 2 

SS 



(nir) 2 yr(2h) 



dx dy 



2tt 



(3.36) 
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which in turn corresponds to the change in coordinates 



x \ ,^(l + ^)Uf r ^]. (3.37) 



vj \ v J V v 

If we consider the theory on a strip of length R then this change in coordinates is an effective 
increase in the length by a factor r n (A). To the order in A to which we are working, this is 
equivalent to a reduction in the strip width by the same factor and consequently a rescaling 
of the eigenvalues of the Hamiltonian by the same factor, r n (A), exactly in accordance with 
equation ( p. 32 ). 



4 Mathematical Arguments 

It is also possible to derive the coupling constant renormalisation and energy rescaling by 
examining the lowest three eigenvalues of the perturbed, truncated Hamiltonian directly. 
This gives exact expressions in terms of integrated four-point functions of the conformal 
field theory. We can then find the large-n behaviour of these eigenvalues by a saddle-point 
approximation. It turns out that these are given in terms of the crossed four-point functions 
and the leading large-n behaviour is given by the leading terms in the crossed four-point 
functions, which are in turn given by the operator product expansion coefficients. In this way 
we see that the "physical" arguments of the previous section are indeed correct. Furthermore, 
we can show that the corrections we have found are the leading order corrections for all the 
eigenvalues of the perturbed Hamiltonian, not just the lowest three. 

We start from the expression (|2.6|) for the dimensionless TCSA Hamiltonian: 



U n = (L - ^) + ~XP n <p(l)P n . (4.1) 

We have introduced A = \ir~ y for notational convenience. 

For simplicity we consider a case where the Hilbert space is a single highest weight rep- 
resentation of the Virasoro algebra with conformal weight H, and that the two lowest level 
states are 

h/>>, £-ih/>>. (4.2) 

We will assume that the the energy-rescaling and coupling constant renormalisation take the 



forms ( 3.13 ) and 

g(x) = l + bTr~ y x + 0(x 2 ,n~ 2y ) , r{x) = 1 + a^ 2y x 2 + 0(x 3 , n~ 3y ) . (4.3) 

We also assume that the first energy gap takes the exact form 

Ai(A) = £i(A) - E (X) = 1 + a\ + /3A 2 + 0(A 3 ) , (4.4) 

so that the TCSA approximation at truncation level n is given by 

A?(A) = r(An^)A 1 (A 5 (An^)) 

= l + a\ + \ 2 (P + an- 2y + ban- y ) + 0(\ 3 ,n- 3y ) . (4.5) 
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This means that we can find the coefficients a and b appearing in the rescaling and renor- 
malisation functions (|4.3| ) from the n dependence in the second order term in the first energy 
gap. From the physical arguments of the previous section, we expect 

By standard perturbation theory, we find the following expressions for the ground state 
and first excited state energies, to second order: 



E 



H + XC' - (AC") 2 
(H + 1) + \C fl4 
+(AC") 2 



1 dz 



(h-l)h 
2H 



-h 2 




2H 


fM 



h 2 /(2H) (2H + h{h-l)f 



{2H) 2 z 1 ~y 



Here C is the boundary coupling constant in the operator product expansion 

+ ... , 



(*m = Si 



and the function F° is the four point chiral block 



(4.7) 



F\z)= ^. 



4> <P 



1 



l (i + --0, 



(4.9) 



that appears in the boundary two-point function on the strip. Since we have chosen boundary 
conditions for which the Hilbert space of the model on the strip has only a single representation 
with highest weight \tp), the boundary two point function is given as the product of the 
structure constants and the single chiral block with the representation t/> in the intermediate 
channel: 

<<Kl)<K*)) S trip = <V#(1)<K*M = (C') 2 F°(Z) . (4.10) 

F 1 is the function that appears in the two point function of the field <p(z) in the first excited 
state, 



F 1 (z) = 

V ' 2H 



\LiL- 

z(z — 



(C'YF^z) , 



l)— + 2zh-l 
dz 



l)— + 2h 

dz 



+ 2H\F 



(4.11) 



If h, the conformal dimension of the perturbing field, is greater than or equal to 1/2 then the 
integrals in the expressions (fOj) diverge, but their difference, the energy gap Ai, is finite. 
We can read off the coefficients a and (3 that appear in (|4.4j ) as 



(X = — — — C , 



2H 



(4.12) 



2H 



F 



1 h 2 /(2H) (2H + h{h-l)f 



(2H) 2 z 1 ~y 



F -r 

z h 



(4.13) 
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These are the exact coefficients in the perturbative expansions of the energy gap; we are 
interested in the n-dependence in the TCSA approximation ( |4.5| ). The n-dependence comes 
from repiacing the functions F° and F 1 in ( 4.13j ) by the TCSA truncations. We see that a 
is independent of n but j3 does depend on n; we shall denote the TCSA truncation by (3(n) 
and we expect from (f4.5[) that 



P(n) =p + an~ 2y + ban~ y + 0(n~ Zy ) . 



(4.14) 



As before, in section |3.4| , the TCSA truncations of the functions F° and F 1 come from 
restricting their Taylor expansions to include only modes up to power z n . It is easier 
to analyse the change in the functions as the TCSA level is increased than to find the n- 
dependence directly, so we consider the change in j3 given by the taking just the coefficient of 



in the integrand of ( |4. 13| ) 



Pin) 



(C 



>Y2 



dz 



2mz n + l 



(4.15) 



The contour in ( 4.15| ) is a small circle around the origin. The function z h (F 1 - F°) is single 
valued around the origin but has a cut from z = 1 to z = oo, and can be expanded in 
increasing powers of (1 — z) . We use the result that 



dz 



2iriz 



n 



a-l 



T(n + a) 
r(a)r(n + 1) ~ r(a) 



(1 + 0(l/n) ) 



(4.16) 



to see that the dependence of the function ( |4. 15| ) on n is determined by the expansion of the 
integrand in powers of (1 — z). The expansion of F° in (1 — z) is determined by taking the 
alternative conformal block expansion of the boundary two point function: 



(C') 2 F°(z) 



id + CC' 

1p 7p 



1> 



+ 



(1 - z)- 2h {1 + 0(1 - z) 2 } + CC'(l-z)' h {l + 0(l-z)} + ... 



(4.17) 



The further terms in ( [4.17 ) correspond to further fields in the operator product expansion of 
4>(x) with itself. We assume that these are less singular than the two shown, as is the case for 
perturbations by the unitary minimal model field 4>i3. Calculating the integral (|4.15| ) we get 



dp _ 2(h-i) 2v _, ,KM , 

dn~ T{2h) H +LL 2HT(h) n + '-' ' 



(4.18) 



where it is worth noting that the leading term which one would expect to have dependence 



n 



2h-2 



n 2y has zero coefficient and the term in n 2y 1 is the subdeading term. We can 
now integrate ( [4.18 ) to find 



Pip) 



/3(oo) + 



T(2h) 



n 



-2y + cc > 



h 



2HyY(h) 



n y + 



(4.19) 
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This exactly agrees with the expected form (4.5), so that the rigorous mathematical derivation 



of the leading n-dependence of the TCSA approximation to the first energy gap agrees with 
the heuristic derivation of the coupling constant renormalisation and energy rescaling from 
physical arguments. 

It is important to note that these effects are independent of the representation ip chosen, 
and so, in particular, they should apply to all of the flows considered by Feverati et al. 

It is in fact possible, by considering the effect of replacing the state L-i\ip) by suitable 
states at arbitrary levels, to extend this calculation to cover all energy gaps, not just the first 
energy gap. This means that if A" (A) is the TCSA approximation to the i-th energy gap, 
then one can prove that 

A?(A) = r(An-^)A,(A 5 (An-^)) + 0(A 3 ) , (4.20) 
where the functions r and g are those given in ([4.3[) with coefficients (|4.6[). 



5 Checks in the tri-critical Ising model 

As mentioned before, we shall perform most of our checks in the tri-critical Ising model 
because the excited state spectra of the strip with a perturbed boundary condition has been 
studied in detail using the TBA approach by Feverati et al. [0, § § allowing us to compare 
the TCSA results with the 'exact' TBA spectrum. The tri-critical Ising model is a unitary 
conformal field theory with central charge 7/10 and the Virasoro algebra has six unitary 
highest weight representations listed in table |]. The model has six fundamental, or "Cardy" , 
conformally invariant boundary conditions |l^, |ll]] listed in table |l[ These can be labelled 
either by representation of the Virasoro algebra or by the allowed values of the boundary 
spins in its realisation as a spin-1 Ising model. The boundary flows were first given by Affleck 
in H and are shown in figure [3|. 



Virasoro label 


(11) 


(21) 


(31) 


(12) 


(13) 


(22) 


Conformal weight 





7 

16 


3 
2 


l 

10 


3 
5 


3 

80 


Boundary spins 


(") 


(0) 


(+) 


(-0) 


(0+) 


(-0+) = (d) 


Boundary fields 


(11) 


(11), (31) 


(11) 


(11), (13) 


(11), (13) 


(11), (13), (12), (31) 



Table 1: The representations of the Virasoro algebra in the tricritical Ising model and prop- 
erties of the corresponding conformal boundary conditions 

From these boundary conditions and their flows one can easily construct superpositions of 
boundary conditions and further flows using the action of topological defects |J. In the case 
in hand, the topological defects and the boundary conditions are labelled by representations 
of the Virasoro algebra, and both the action of the defects on the boundary conditions and 
the representation content of the strip Hilbert space are given in terms of the fusion of 
representations. If we denote the representations of the Virasoro algebra by a, b etc, the 
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Figure 3: The space of boundary flows in the tricritical Ising model 

defects by D a and the boundary conditions by B a then the defects act as 

D a -B b = B a * b = ^ N ab c B c , (5.1) 

c 

where N ab c are the Verlinde fusion numbers. Likewise, the representation content of a strip 
with conformal boundary conditions a and b, which we denote "H( a ,b) is 

«(o,6) = ®cN ab c L c , (5.2) 

where L a is the highest weight representation with label a. The commutativity of the fusion 
algebra means that the strip with boundary conditions (a, b*c) has the same Hilbert space as 
the strip with boundary conditions (b* a, c), which we can interpret as the statement that the 
action of a defect of type b on either of the two boundary conditions a or c leaves the Hilbert 
space unchanged. This means that the six separate flows considered by Feverati et al. can 
now be group into three pairs which can be found as the spectrum of the strip with the basic 
set of flows fl2.7| ) on one edge and one of the three 'fixed' conformal boundary conditions (rl) 
on the other. By considering the second boundary conditions as the action of the defect -D( r i) 
on the (11) boundary condition, they have alternative interpretations as the spectra of strips 
with perturbed boundary conditions of type (12), (13) and (22) coupled with an undeformed 
(11) boundary condition on the second edge as we show in table ||. 

5.1 The coupling constant renormalisation 

In the tri-critical Ising model we first consider the renormalisation of the coupling A of the 
perturbing field 4>3/5 which generates the flows 

(11) & (12) ±3 (21) . (5.3) 

We recall that there are two predictions from this calculation - firstly the general prediction 
that the coupling renormalisation has a particular scaling form, namely that the function g n 
in equation |2.9| has the form 

9n(x) = g{xn~ y ) (5.4) 
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Spectral flow 


Strip configuration 1 


Strip configuration 2 


(11) ^ (12) _). (21) 
(31) ^ (13) (21) 
(21) (22) -> (11)+(31) 


(11; 11) <- (12; 11) -> (21; 11) 
(31; 11) <- (13; 11) -> (21; 11) 
(21; 11) <- (22; 11) -> (11+31; 11) 


(11; 11) <- (12; 11) -> (21; 11) 
(11; 31) <- (12; 31) -> (21; 31) 
(11; 21) <- (12; 21) -> (21; 21) 



Table 2: The possible interpretations of the flows considered by Feverati et al. Configuration 
one has the boundary condition (11) on one side of the strip; configuration two has a fixed 
boundary condition of type (rl) on one side and the basic flow (11) <— (12) — > (21) on the 
other. 



and secondly a particular prediction for the 1-loop behaviour, 



Cx 



The values for the tri-critical Ising model are 

C = - )"' y = 0.544542... , h = 3/5 , y = 2/5 . (5.6) 

For the TCSA approach, we need to specify two boundary conditions. On one we of course 
take the perturbed boundary conditions ( |5.3|) ; on the other we can take any of the "fixed" 
boundary conditions (rl) and we will obtain one of the cases investigated by Feverati et al. 
For simplicity, we start with the boundary condition (11) on the other. For the flow with 
negative A, that is for (11) <— (12), we plot the functions log(<? n (— xn v )) against log(— x) 
for several values of n together with the 1-loop prediction for \og{g{—x)) in equation (|5.5|). 



These are shown in figure 4(a) We see that the functions \og{g n {— xn y )) are indeed almost 
identical for various values of n and in good agreement with the prediction, confirming both 
the prediction of the scaling form and the approximate numerical expression for this scaling 
function. For positive values of A, we plot \og{g n {xn y )) against log(x) for the same values of 
n, again together with the 1-loop prediction for \og(g(x)) in equation ( |5.5| ) in figure 1(b) with 



the same results for small values of A. 

The scaling form does clearly break down for larger values of A positive, where the func- 
tions g n (xn y ) have n-dependent maxima at n-dependent values of x. This corresponds phys- 
ically to the boundary condition approaching close to the (21) fixed point but then moving 
away in the direction of the (13) fixed point. Presumably the flow in this region is governed 
by the critical exponents around the (21) fixed point and different arguments are required to 
analyse this behaviour, but we can make some suggestions about the energy-level-dependence 
of the position of the fixed point, which we do later in this section. 

Finally we can check whether the coupling constant renormalisation and energy re-scalings 
depend on the second boundary condition on the strip. In figure 5(a) we present the coupling 



constant renormalisation as calculated from the three different choices of second boundary 
condition, (11), (21) and (31). We see that, modulo numerical inaccuracies, the three different 
TCSA strip configurations give the same coupling constant renormalisation, as would be 
expected on physical grounds, and that this is in agreement wit the 1-loop calculation. 
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0.30 



-4.0 -3.5 - -2.5 -2.0 -1.5 -1.0 PD 




(a) The scaling functions log g n (— xn v ) and the pre- (b) The scaling functions log g„(xn v ) and the predic- 
diction for \og(g(— x)) plotted vs. log(— x) i.e. for tion for \og(g(x)) plotted vs. log(:c) for positive cou- 
negative coupling constant pling constant 



Figure 4: The numerical coupling constant renormalisation found for the tri-critical Ising 
model with (11) boundary condition on the other edge, for truncation levels 8 (o) 15 (•) and 
22 (□) together with the 1-loop predictions ( 3.14| ) shown as a solid line. 



5.2 The energy rescaling 



We can also check the prediction 3.32 for the form of the energy rescaling in the tri-critical 
Ising model. This is shown in figure |5(b) where log(r(— x)) is plotted together with the 
numerical estimates for log(r n (— xn y )) for various different truncation levels and for different 
choices of the non-perturbed boundary condition. This confirms both the scaling form of v n 
and also the numerical coefficient calculated in this section. 



6 Checks in the Yang-Lee model 

The Lee- Yang model two conformal boundary conditions which are connected by an integrable 



boundary flow which has been studied in great detail [14|. The model is nonunitary and the 
boundary field generating the boundary flow has conformal weight —1/5, which means that 
the perturbation is UV finite but IR divergent. Since h = —1/5, y = 1 — h = 6/5 and 
consequently the leading corrections from truncation errors, which are power series in n~ y , 
decay very rapidly and the TCSA quickly becomes very accurate. The truncation does still 
have an effect, of course, and the accuracy of the numerical method allows us to check the 
predictions for the coefficients of A 2 in the first and second energy gaps. The results are shown 
in ^| where we show the coefficients extracted from the TCSA method and also the prediction 
based on the two leading terms in n~ v . We see very good agreement so that the leading terms 
do indeed give the correct behaviour. We have also checked this in the tri-critical Ising model 
where we also get good agreement, but the corrections from the terms in n~ 3y are larger. 
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0.10 - 



-0.10 - 




0.06 




(a) The TCSA coupling constant renormalisation (b) The energy rescaling function log r„(— xn v ) and 

log(<7n(— xn y )) plotted against log(— x) for the tri- the prediction for log(r(— x)) plotted vs. log(— x) for 

critical Ising model and choices of second boundary negative coupling for the flow (11; 11) <— (12; 11) for 

condition (11) = (-), (o, level 22), (21) = (0) (•, truncation levels 8 (o) 15 (•) and 22 (□). 
level 14) and (31) = (+) (□, level 18) together with 
the 1-loop prediction (3.14) shown as a solid line. 



Figure 5: The coupling constant and energy rescaling functions 



We have used the following expression for the second energy level in the Lee- Yang model, 
which is valid for any system where the representation ip of conformal weight H is degenerate 
at level 2 and has only a single state at that level which we have chosen to be L_2\ip)'- 



En 



= + 2) + AC" ( H 
+ (AC') 2 

z(z 2 



4h(h-l) 



l 



2h 2 h 2 (2h-\) 2 
AT + 2HM 



^1 

2" 



l)— + {'iz 2 + l)h-2 



Ah 2 h 2 (2h-l) 2 {Ah{h-l)+M) 2 

{M) 2 z l -y 
+ m) f° , 



Mz 3 -y (2HN)z 2 -y 

(z-l/z)— + h(3 + l/z 2 ) 





1> 




1 z 



4> 



+ 4H , 



where for the lee- Yang model h = H = — 1/5 and 



r(i/5)r(2/5) \ 
r(-i/5)r(4/5); 



1/2 



(6.1) 



(6.2) 
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Figure 6: The coefficients of A 2 in the first and second gap for the Yang-Lee model. The dots 
are from TCSA data, the solid curve the prediction using the two leading terms in 1/n. The 
straight line is the exact value. 



7 Fixed points and the coupling constant renormalisation. 

One of the principal effects of the coupling constant renormalisation in the TCSA approach is 
the possibility to move the IR fixed point of the finite-size scaling flow from infinite coupling 



constant to finite coupling constant. From equations fl2.9|) and ( 3.14 ), we see that for C 
positive, the IR fixed point at Aqo = oo is brought to the finite value 



a: 



c 



{irn) 1 



(7.1) 



However, just as the coupling to the identity operator depends on the level of the state, so 
the coupling constant renormalisation depends on the level of the state, so that the one-loop 
effective coupling constant renormalisation of a state at level E is not given by (3.13), but 
rather by 

X* n {E)=A(n){n-E)y , (7.2) 

where the one- loop prediction for A(n) is A{n) = yT(h)TT y jC . The exact position of the fixed 
point is not determined very accurately by the 1-loop calculation, (indeed we shall see in 
section |8.2| that the position of the fixed point in the tri-critical Ising model is approximately 
A* = By/n) but we might conjecture is that the dependence on energy level still has this 



form. In figure 7(a) we show the full TCSA spectrum for the tri-critical Ising strip with 
boundary conditions (12) and (11). The level crossings seen for positive A occur at the IR 
fixed point, and their position depends on the excitation level. We have also shown a fit to 
these positions of the form Afp {n,E) = A(n)(n — E) y . This at least appears to capture the 
qualitative behaviour quite well. 

A second prediction is that the position of this fixed point is independent of the unper- 
turbed boundary condition. In figure 7(b)| we show the TCSA spectra of the three strips 
with one boundary being the perturbation (12) —> (21) and the other being one of the three 
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fixed boundary conditions (rl). As can be seen, the three spectra agree very well indeed 
on both the location of the fixed points, indicated by the multiple line degeneracies, on the 
energy-dependence of these fixed points and on the rescaling of the Hamiltonian, confirming 
the prediction of our analysis that these effects are independent of the unperturbed boundary 
condition (to leading order). 




0.5 1.0 1.5 2.0 

(a) The full TCSA spectrum showing the variation 
of the fixed point with energy level; also shown is 
the best fit to these positions in the scaling form 
Ad-it = A(n — E) v , in the case n = 14. 
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5 



4 



3 



2 
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-4 -3-2-10 1 2 

(b) A simultaneous plot of the TCSA spectrum of 
the basic flow (12) — > (21) on one edge and the fixed 
boundary condition (11) (red), (21) (blue) and (31) 
green on the other edge for truncation level 14. 



Figure 7: Energy-dependence and boundary-condition independence of the position of the 
fixed point in the flow (12) — > (21). 



8 Flows generated by irrelevant fields 

When we consider a flow with a UV and an IR fixed point, we can attempt to describe it in 
terms of perturbations of either fixed point. A perturbation of an IR fixed point will be by an 
irrelevant field which is a non-renormalisable perturbation. This means that the perturbation 
theory will have divergences that require regulation and that a possibly infinite number of 
counterterms will be required to allow one to remove the regulator. For this reason, they 



have only infrequently been studied in the literature [18, 19, pQ] . In TCSA, however, there is 
no need to introduce counterterms, one can simply compute the spectrum of the truncated 
Hamiltonian. 

We have investigated this in the case of the tri-critical Ising model. There are three IR 
fixed points in the basic sequence (fO|). In the case of the (11) and (31) boundary conditions, 
the fields on the boundary comprise the vacuum representation of the Virasoro algebra so 
it is expected that the flow will be generated by the field T(x), the boundary stress-energy 
tensor, as in [p0|]; for the (21) boundary condition, the boundary fields form the (11) and (31) 
representations and it is to be expected that the flow is generated by the field </>(3i) of weight 
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3/2, mirroring the case in the bulk flows considered in |T{|]. We shall present some results for 
this last case, and for the flow (21) — > (13), in particular. 



8.1 Comparison with exact results 



Firstly, we find that the TCSA spectrum is in very good agreement with the "exact" TBA 
spectrum; we show the results for the spectrum in figure |ll[ Secondly, we can compare the 
TCSA couplings A^ 21 ^ at truncation level n with the effective coupling A^ 3 ^ for the description 



of the system as a perturbation of the (13) by the field 
expect a relation of the form 



(13V 



On dimensional grounds, we 



a£ 3) = (A( 2i >r 



4/5 </(Ai 21) 



n 



1/2 



and this indeed what we find. In figure S(a) we plot the effective couplings against the 



TCSA coupling for various truncation levels, and in figure 3(b) we plot the scaling function 




Figure 8: Truncation level dependence and scaling functions for the irrelevant coupling A^ 21 ) 
at truncation levels 8 (red, dotted), 13 (green, dashed) and 18 (blue, solid). 



The scaling form Q8.1| ) does, however, have the major consequence that the TCSA coupling 
constant does not converge to a fixed value as the truncation level is increased. For an 
irrelevant perturbation, h > 1 and so y < 0. If the function g(x) in ( 3.13| ) is approximately 



equal to one for a range |x| < A, then the TCSA scaling function g n (A) is approximately equal 
to one for |A| < An" y . For y positive, this range grows with increasing n; for an irrelevant 
perturbation it shrinks. In this case, we see that 

A£ 3) (A( 21 >) 4 / 5 ^0.06 for A( 21 )<^. (8.2) 



8.2 Comparison with TCSA parameters 

While (|8.l|) and (|8.2|) relate the TCSA parameter An 21 '' to the "exact" parameter A^ , it is also 
interesting to compare the two TCSA parameters. In this case we find the problem that the 
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ranges in which the TCSA parameters are good approximations to the "exact" parameters do 
not overlap. The IR parameter is good for A( 21 > < 0.05 n- 1 / 2 , which equates to A^ 13 ) > .7n 2 / 5 , 
while we can see from figure 4(a) that the UV parameter is good for A^ 13 ) < .5n 2 / 5 (where we 
have used the Z2 symmetry to relate A*- 13 ) to — A*- 12 ).) As a consequence, there are no values 
of the TCSA coupling constants for which (8.2) holds. Instead we find the interesting plot 



p(a) ; as n increases, the relation between the two TCSA couplings changes. For either of the 



couplings small, there is an approximate linear relationship. We can elucidate this by noting 



that the TCSA coupling relation itself takes a scaling form as we show in figure 9(b), where 
we see that, to a good approximation, 



(8.3) 



where 



F(x) ~ 0.5 — 6.x for x small, F(x) ~ 0.3 — 1.3x for F(x) small 



U) 



The linear relation is to be expected, as the perturbing fields at one end (UV or IR) will flow 
into the perturbing field at the other (IR or UV, respectively) and so the coupling at one end 
will end up being linearly related to the coupling at the other. Quite why this relation takes 
the form ( |S.3D is unclear at the moment. 




9 TCSA flows beyond the fixed point 



As has been remarked on before Jl2| , 2C], the TCSA spectra have the remarkable property 
that they extend beyond the IR fixed points and can encompass several different perturbative 
flows. These appear to follow the sequence of flows found first by Lesage et al in [15] and 
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-0.15 -0.10 -0.05 



0.05 0.10 0.15 



Figure 10: The scaled gaps for the perturbation of the (21) boundary condition on a strip 
with boundary conditions (21; 11) by A^ 21 V(3i)) plotted against A^ 21 -* at truncation level 14. 
The four dashed lines are the approximate positions of the fixed points (in order) (11), (12), 
(13) and (31). 



which also appear in Fredenhagen et al |16| and Dorey et al |T7[] and which in the simplest 
form applicable to the tri-critical Ising model is 

(12) — » (11) <- (12) -> (21) <- (13) -»• (31) <-— (13) (9.1) 

where the dashed arrows reflect what is seen in TCSA and are a natural extension of the 
sequence in the papers cited. The same sequence is seen no matter which boundary condition 
we take as the starting point from which we perturb. This applies also to the case of irrelevant, 
non-renormalisable perturbations starting from the 'fixed' (rl) boundary conditions. This 
enables us to put coordinates on the sequence, in the sense that every pair of flows is covered as 
the standard perturbation of a boundary condition, and we can relate the coupling constants 
on successive overlapping pairs: 



(12)+/i( 12 )0 13 (13)+/z( 13 V 



13 



(12) — > (11) <- (12) -> (21) «- (13) -> (31) *- (13) (9.2) 

V V V 

(U)+^T+.. (21) +/ u( 21 )0 31 + .. (31) +/ u( 31 )T+.. 

where the ellipsis in the perturbations of the fixed boundary conditions shows that further 
counterterms may be needed as these perturbations are non-renormalisable. As an example 
of the way the TCSA flows can encompass several fixed points, in figure || we show the scaled 
gaps for the perturbation of the (21) boundary conditions by the irrelevant field <^>(3i) of 



conformal weight 3/2. As can be seen, all five fixed points of (9.1) appear. 

To give support to this picture, we present the TCSA spectra for the flow (13) — > (21) 
obtained from the three separate routes: the relevant perturbation of the (13) boundary 
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Figure 11: The scaled gaps for the flow from the (13) boundary to the (21) boundary condition, 
as found from the relevant perturbation of the (13) boundary condition (•), from the irrelevant 
perturbation of the (21) boundary condition (o), and from the extension of the flow generated 
by the relevant perturbation of the (12) boundary condition (□) at truncation level 22. They 
are plotted against £, the boundary parameter in Pearce et al. 



condition, the irrelevant perturbation of the (21) boundary condition, and the continuation 
of the flow (12) — > (21) given by the relevant perturbation of the (12) boundary condition. 



These are all shown in figure 11 



The TBA spectra are given as functions of the boundary reflection parameter 6 b in 
£ of Feverati et al 



23] 



Iflor 

we can find the exact 



J] (note £ = 9b)- Using the results in J22| 
expression for 6b in terms of the perturbative coupling constant of the relevant pertur 
bation of the (13) boundary condition. We can also find the approximate numerical result for 
the boundary parameter in terms of the TCSA coupling A^ 21 \ which is only valid for small 
coupling. These are summarised below 



Boundary 


Ob 


(13) 


(f log a- log d) + |logA( 13 ) 




= 1.96102... + flog A^ 13 ) 


(21) 


-5. - 2.1ogA( 21 ) + . . . 



(9.3) 



The constants in the first relation are 



§7rsin(§7r) 



r(i)3r(!)sin(^)sin(f7r) 



(9.4) 



2 The equations in [jLEiL for the ^-function are known not to be correct, but they do describe the changes 
under a purely boundary perturbation. See j^J for the correct equations when there is a simultaneous bulk 
perturbation. 
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Finally, we have also found an approximate numerical relation between the couplings A^ 21 -* 
and A^ 12 ) for this region, which is 

A (12) 

-2=- - Z.jnX^ ~ 0.5 , (9.5) 

which is in agreement with ( |8.4| ) to the accuracy that is obtainable for the flow from the (12) 
boundary condition beyond the (21) fixed point. 



10 Conclusions 

We have made further progress in understanding the errors in the truncated conformal space 
approach to perturbed boundary conformal field theory on a strip. The principal quantity 
calculated with TCSA is the spectrum as a function of the coupling constant. The errors are 
of two principal sorts - a renormalisation of the coupling constant and a change in the energy 
scale. We have found perturbative expressions for these from physical arguments and showed 
they are correct using an analysis of the perturbative spectrum. 

One important aspect of these predictions is that they are independent of the second 
boundary condition on the strip and we have verified this by considering different choices for 
the unperturbed boundary condition and finding the same results for each choice. 

Furthermore we have investigated the behaviour of the TCSA spectrum for large coupling 
constant. We have reported before that this appears to show a sequence of RG flows in the 



same pattern as found by Lesage et al [15|, and we know show that this is quantitatively 
correct as well, in the case of the tri-critical Ising model. 

It would be good to find some way to describe the full RG flows such as in figure || in 
terms of some beta-function which has the sequence of fixed points as its zeroes, but at the 
moment we are unsure how to do that. In the lattice model, the perturbation parameter is a 



boundary magnetic field and the sequence of fixed points (9.1) (excluding the points joined 



by the dashed lines) can be found simply by varying this field |2^]; of course in the quantum 
model this is not simply the case - the natural sequence of flows splits up into overlapping 
pairs of flows in both the TBA and the TCSA descriptions. It may be helpful to use the 
coordinates on the full moduli space given by the values of the ^-function and the excited 
(^-function which corresponds to the overlap with the bulk spin field. This is something we 
hope to return to shortly, 

Finally the ideas on the perturbative treatment of the corrections to the TCSA presented 



here can be easily adapted to the case of bulk flows which we plan to address in [25] 
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